﻿ Vladimir Kotelnikov John Wozencraft Irwin Jacobs Simon Haykin 1 Detecția mai multor semnale din zgomot Abordarea detecției de până acum este cea a unui operator de RADAR sau SONAR În ipoteza H0 se presupunea că doar zgomotul este prezent, nu și semnalul util (cel căutat) În ipoteza H 1 se presupunea că este prezent semnalul util, dar afectat de zgomot În cazul sistemelor de telecomunicații abordarea detecției este diferită În mod uzual se transmite unul din cele M semnale, corespunzătoare celor M simboluri cu care se fac transmisiile La recepție nu se pune problema să decidem dacă este prezent un semnal util sau numai zgomotul Trebuie să detectăm care dintre cele M semnale este prezent, evident, însoțit de zgomot aditiv Deși ceeace descriem este o problemă de clasificare, o vom numi tot detecție Vom discuta cazul detecției a M=2 semnale și vom trece apoi, prin generalizare, la cazul M>2 Detecția unui semnal dintre două semnale posibile Definim problema de detecție, sau, dacă doriți, de testare a ipotezelor +H 00: [ ] [ ] [ ]; n=0, 1, , N-1xn s n wn= H 11: [ ] [ ] [ ]; n=0, 1, , N-1xn s n wn=+ 2 , σw0I∼ () uN Cele două semnale utile se consideră a fi deterministe, și pe deplin cunoscute, iar zgomotul este alb, gaussian În telecomunicații erorile de orice tip sunt la fel de indezirabile, motiv pentru care vom face detecția recurgând la criteriul bayesian al probabilității de eroare minime Dacă cele două simboluri utilizate în transmisie sunt echiprobabile, tot așaeste și apariția celor două semnale utile, așa că HH {}{} 0112PP== Detecția bayesiană se face conform valorii raportului de plauzibilitate condiționată, adică conform relației (1) 1H H {} {} 10pPxH 1=> {} {} 01pPxHH ceeace înseamnă că detectorul este de tip ML Cu datele x recepționate, se calculează valorile ()() {} 01;ppxxHH Decizia se ia în funcție de care dintre aceste două valori este mai mare Repartițiile datelor în cele două ipoteze se deduc ușor șisunt 2 , ; 0, 1iσ=xsI∼N () iu adică 1N− 211⎧⎫ xH −= () () ⎨⎬ ∑ 22exp [ ] [ ] ; 0, 1iNipxnsni=− 2 022nσ = ⎩⎭ πσ () 2 Maximizarea expresiei ; 0, 1pi=xH () i este echivalentă cu minimizarea exponentului, mai precis a sumei 1N− 22 Dxnsni=− =[ ] [ ] ; 0, 1 (16) () ii∑ 0n= Cu datele x[n] “măsurate”, adică recepționate, și cu cele două semnale utile păstrate în copie (replică), se calculează, aplicând relația (16), valorile 22 {} 01;DD și se alege acea ipoteză care corespunde valorii minime Reamintim că, dacă u este un vector în spațiul real N dimensional, pătratul normei sale este 1N− 2 2T uuu ∑ ku== 0k= În spațiile cu dimensiune N cel mult 3, norma vectorului este egală cu lungimea sa (euclidiană) Vom păstra pentru normă denumirea de “lungime”, deși ea poate fi improprie Distanța dintre punctele definite de doi vectori, u și v,se calculează cu relația 1N− 2T 2 ,duv=− −= −uv uv uv ()()() () kk∑ 0k= Vom privi vectorii xw s s , 01,, ca fiind din spațiul real, N dimensional După cum se vede dintr-un exemplu prezentat selectarea ipotezei H1 în figură, pentru un spațiu bidimensional, D1 ⋅+>⋅xx xxHH () ()() () {}{} 100 011 12 12PT T PT T= Se introduce o nouă statistică, T(x), ca diferență a celor două statistici utilizate 1N− 1 = −−−xxxεε () () ()()() ∑ TTT xnsnsn=−1010 10[ ] [ ] [ ] (17) 02n= Datele x au, în ambele ipoteze, o repartiție gaussiană 2 , ; 0, 1iσ=xsI∼N () iu iar statistica T(x) are o dependență liniară de aceste date, așa că, în ambele ipoteze șiT este gaussiană Este suficient, deci, să determinăm mediile și dispersiile ei, în ambele ipoteze, pentru a scrie expresiile plauzibilităților Pentru calculul mediilor avem 5 1N− 1⎧⎫ −− HHεε ()() {} ⎨⎬ ∑ 10010 0[] [] []ET E xn s n s n=− 02n= ⎩⎭ 1N− 1 εε − H ()() {} ∑ = [ ] [ ] [ ]Exn sn sn−−1001 0 02n= 1N− 1 −−− εε {}( )() ∑ [ ] [ ] [ ] [ ]Es n wn sn s n=+10010 02n= 1N− 1 −− εε ()() ∑ [ ] [ ] [ ]sn sn sn=−1001 0 02n= 0ε 11NN−− 21 () −−εε () ∑∑ [ ] [ ] [ ]snsn sn=−1001 0 002nn== și rezultă media condiționată a statisticii T 1N− 1 +Hεε () {} ∑ 10010[] []ET s ns n=− 02n= Această formă mai poate fi rescrisă astfel 1N− 1 H +εε () {} ∑ 10010[] []ET sns n=− 02n= − ⎧⎫ 211 121NN N−− − + ()() ⎨⎬ ∑∑ ∑ 1100 [ ] 2 [ ] [ ] [ ]sn snsn sn=− = 00 02nn n== ⎩⎭ 1N− 21 − () ∑ 10 [ ] [ ]sn sn=− 02n= 21 − ss =−10 2 În mod asemănător se arată că 21 ssH {} 110ET =− 2 Dispersiile se calculează cu 1N− 1⎧⎫ −− εε HH ()() {} ⎨⎬ ∑ 100100[] [] []Disp T Disp x n s n s n=− 02n= ⎩⎭ 11NN−− 22 2 H −= − ()() {} ∑∑ = [ ] [ ] [ ] [ ] [ ]Disp x n s n s n s n s nσ01 0 1 0 00nn== 2 2 ss =− σ10 Dispersia în cealaltă ipoteză este aceeași 2 2 =−ssH {} 110Disp Tσ 6 Am ajuns în etapa în care putem descrie repartițiile statisticii T, în cele două ipoteze 221⎧ 2⎛⎞ H −− − ss ssN 10 10 0, σ , în ipoteza ⎜⎟ ⎪ 2⎪⎝ ⎠ T⎨ ∼ 221⎛⎞ 2⎪ H ss ssN +− − 10 10 1, σ , în ipoteza ⎜⎟ 2⎪ ⎝⎠ ⎩ În figură se arată repartițiile statisticii T, în cele două ipoteze Densitățile de probabilitate diferă doar prin medie, dispersia fiind aceeași Analizând simetria figurii, se remarcă ușor că cele două tipuri de erori ce apar în cursul luării deciziilor, sunt egale, adică = Dar 22⎛⎞ ⎛⎞ ss ss − 10 100121+− ⎜⎟ H >= = {} 020PT Q Q⎜⎟ 2 22σ ⎜⎟ ⎜⎟ ss − 10σ ⎝⎠ ⎝⎠ și 2 ⎛⎞ ss 101− ⋅ xx xxHH () ()() () {}{} 100 011 12 12PT T PT T=> ⎧⎫⎧⎫ ⎪⎪⎪⎪ >⋅+− 01 2 ⎛⎞ ss 101− ⎜⎟ (18) Q=2 2σ ⎜⎟ ⎝⎠ Este evident că, pentru a reduce acest minim al erorii, trebuie crescută distanța dintre cele două semnale cu care se face transmisia, dar respectând constrângerile energetice ce apar în telecomunicații în general și în cele mobile în special 7 Energia medie (statistic) utilizată într-o transmisie binară, cum este cazul analizat, se calculează cu relația εε 01+ =HHε εε {}{} 0101PP+= 2 Pătratul normei diferenței vectorilor de semnal se poate pune și sub forma 2T TT T − −= + − ss ss ss ssss ss ()() 10 10 10 1100 102−= T ⎡⎤ 10TTss − ss ss ⎢⎥ () 11 0 0 1TT=+ ss ss () ⎢⎥ 11 0 012+ ⎣⎦ T ss 10⎛⎞ εε − () 10 1=+ ⎜⎟ ε ⎝⎠ Se definește coeficientul de corelare al celor două semnale cu care se transmite, cu relația TT 10 10ss ss ; 1 1==−≤≤ρρ TTss ss s sε () 11 0 02+ Dacă cei doi vectori de semnal sunt ortogonali, coeficientul lor de corelare este nul; dacă vectorii de semnal sunt coliniari, dar au sensuri opuse, valoarea coeficientului de corelare este -1 Relația anterioară devine 2 −= −ssε () 1021sρ Substituim acest rezultat în forma (18) a probabilității de eroare, șiobținem pentru ea ⎛⎞ 1ρε () s− = 2e (19)PQ⎜⎟ ⎜⎟ 2σ ⎝⎠ În concluzie, pentru a obține cea mai redusă valoare a probabilității de eroare, la o energie medie dată și o putere a zgomotului dată, trebuie să semnalizăm cu semnale coliniare, dar de sensuri opuse Astfel de perechi de semnale se numesc semnale “antipodale” Exemplu Sistem de transmisie binară, cu modulaţie coerentă de fază, PSK Într-un sistem PSK binar se transmite unul dintre semnalele antipodale, având coeficientul de intercorelaţie −1 =∈ ()() 000[ ] cos 2 ; 0, 0 5sn A fn fπ π=+=− ()() snA fn A fnππ10 0[ ] cos 2 cos 2 ; =− =− ⇒ 10[] [] ssn snρ 1; Pentru precizarea ideilor, vom considera o frecvenţă digitală de 0 25 cicli/eşantion, pentru care semnalele utilizate în transmisie devin π ⎛⎞ sn A n n==0[ ] cos ; 0, 1 ⎜⎟ 2⎝⎠ π ⎛⎞ = 1[ ] cos ; 0, 1sn A n n=− ⎜⎟ 2⎝⎠ =− =− ⇒ 10[] [] ssn snρ 1; 8 Cei doi vectori de semnal sunt acum AA− ⎡⎤⎡⎤ ss 01; ;== ⎢⎥⎢⎥ 00 ⎣⎦⎣⎦ şi sunt reprezentaţi în figură Se vede că sunt vectori coliniari, dar de sensuri opuse Vom calcula în mod direct coeficientul de intercorelaţie, deşi cunoaştem valoarea sa A⎡⎤ [] ⎢⎥ 220TA− ss ⎣⎦ 10022A− ====− 221TTsρ ss ss ++ ⎡⎤ ⎡⎤ 11 0 0AAAA− 00AA−+ [] [] ⎢⎥ ⎢⎥ 00 ⎣⎦ ⎣⎦ Probabilitatea minimă de eroare devine, conform relaţiei (19), ⎛⎞ ε ⎜⎟ 2e ; (20PQ=) σ ⎝⎠ ,===εεεε pentru un simbol binar 10unde este energia emisă Pentru un raport ENR de 15 dB (adică 31 6 în valoare absolută), rezultă o probabilitate medie de eroare de 8− () e31 6 10PQ≅≅ Curba probabilitate de eroare funcţie de ENR se dă în figură, împreună cu curba corespunzătoare sistemului FSK Exemplu Sistem de transmisie binară,*** cu modulaţie coerentă de frecvenţă, FSK Într-un sistem de transmisie binară, cu modulaţie FSK coerentă, se efectuează transmisia cu semnalele ortogonale =∈= ()() 000[ ] cos 2 ; 0, 0 5 , n 0, 1, , -1sn A fn f Nπ =∈= () () 111[ ] cos 2 ; 0, 0 5 , n 0, 1, , -1sn A fn f Nπ 1 ⇒⇒ −== 001Pentru: [ ] [ ] 0 sff snsnρ 0; 2N Cele două semnale au aceeaşi energie, energia emisă pentru un simbol binar, adică ==εεεε 10= şi deci probabilitatea minimă de eroare este, conform relaţiei (19) ⎛⎞ ε ⎜⎟ 2e; (21PQ=) 2σ ⎝⎠ 9 Din (20)şi (21) vom deduce relaţia dintre puterile necesare pentru transmiterea unui simbol binar, cu aceeaşi probabilităţi de eroare, în sistemele PSK şi FSK ⎛⎞⎛ ⎞ ε ε PSK FSK = εε ⎜⎟⎜ ⎟ FSK PSK 22e= 2PQ Q=⇒ 2σσ ⎝⎠⎝ ⎠ ε ε ε FSK PSK PSK =+ 10log 10 log 10 log 2 10 log 3dB=+22 2sau σσ σ Ca să asigure aceeaşi probabilitate de eroare, un sistem FSK are nevoie de o putere emisă, pe simbol binar, dublă faţă de cea necesară unui sistem PSK Curbele probabilităţilor de eroare, ca funcţie de ENR, atât pentru sistemul PSK, cât şi pentru sistemul FSK, sunt arătate în figură Se remarcă diferenţa de 3dB între cele două curbe *** 10 Receptoare M-are Dacă avem de transmis nu două, ci M>2 simboluri, vom folosi pentru transmisie M semnale, câte unul pentru fiecare simbol {} 01 1[], [], , []Msnsn s n− În telecomunicaţii, cele M simboluri au probabilităţile apriorice egale între ele şi deci de valoare 1/M Receptorul (detectorul) va fi deci unul de tip ML Avem regula de decizie Pentru 1 , 0, 1, , -1PiM== {} iH M cu datele x recepţionate, se calculează cele M valori ale plauzibilităţilor condiţionate , 0, 1, , -1piM=xH {} i Dacă , 0, 1, , -1 este ipoteza adevărată, (22)ppiki M=⇒>≠xxHH H {} {} ikk *** Receptorul ML ce funcţionează conform regulii (22) este un receptor optimal El este tot un receptor de distanţă minimă Schema sa este prezentată în figură şi constă din M corelatoare de replică, ce generează M statistici şi dintr-un bloc de selectare a valorii maxime Se calculează statisticile 1N− 1 xε [ ] [ ] ; 0, 1, , -1Txnsni M=− () ∑ iii= 02n= şi se ia decizia în funcţie de care anume dintre valorile calculate este mai mare 11 Determinarea probabilităţii de eroare este dificilă în acest caz, deoarece apare o eroare dacă oricare dintre cele M−1 statistici depăşeşte statistica asociată cu ipoteza adevărată Dacă semnalele cu care se efectuează transmisiile {} 01 1[], [], , []Msn sn s n− sunt ortogonale între ele, atunci statisticile sunt şi ele ortogonale între ele, şi sunt şi necorelate Deoarece datele sunt gaussiene, şi statisticile, care sunt funcţii liniare de date, sunt gaussiene iar necorelarea atrage după sine independenţa (statistică) a statisticilor Pentru a verifica acest fapt calculăm covarianţa a două statistici Mai întâi 1N− 1⎧⎫ [][]ExnsnET− HH ε = {} ∑ lliii⎨⎬ 02n= ⎩⎭ 1N− 1⎧⎫ ε [ ] [ ] [ ]Esnwnsn=+ () ⎨⎬ ∑ iil− 02n= ⎩⎭ ⎧⎫ 11NN−− ⎧⎫ 1⎪⎪ [ ] [ ] [ ] [ ]EsnsnE wnsn=+ ε ⎨⎬⎨⎬ ∑∑ iiil− 002nn== ⎩⎭ ⎪⎪ 0= ⎩⎭ 1N− 1 [ ] [ ]Ewn sn=− {} i∑ iε 2 0n= 1 ; , 0, il M i l=− 1, , 1; − ≠ ε i = 2 Calculăm covarianţa, dacă i şi j diferă de l =− − HH HH ,ET ET T ETCovTT⎡⎤ {} {}{} {} ljlljjliii⎡⎤ ⎣⎦ ⎣⎦ 11NN−− 11 1 1⎧⎫ ⎡⎤⎡ ⎤ −+ H εε ε ε [ ] [ ] [ ] [ ]E xns n xns n=− ⎨⎬ ∑∑ liii j jj+ ⎢⎥⎢ ⎥ 0022 2 2nn== ⎣⎦⎣ ⎦ ⎩⎭ 1N− 1N− ⎧⎫ ⎡⎤ ⎡⎤ []n⎨⎬ [ ] [ ] ] [ ] [ ][Esnwnsn snnsw=+ ()() ∑ ∑ ijll + ⎢⎥ ⎢⎥ 0n= 0n= ⎣⎦ ⎣⎦ ⎩⎭ 11NN−− ⎧⎫ ⎡⎤⎡ ⎤ [ ] [ ] [ ] [ ]E wns n wns n= ∑∑ ij⎨⎬ ⎢⎥⎢ ⎥ 00nn== ⎣⎦⎣ ⎦ ⎩⎭ 11 11NN NN−− −− ⎧⎫ [ ] [ ] [ ] [ ] [ ] [ ] [ ] [ ]Ewnwms ns m E wnwm s ns m== {} ∑∑ ∑∑ ij ij⎨⎬ 00 00nm nm== == ⎩⎭ 11NN−− 2 [ ] [ ] [ ]nmsnsmσδ ∑∑ ij=− 00nm== 1N− 2 [][]snsiji lnljσ=≠≠≠=0; pentru , si ij∑ 0n= 12 Vom mai simplifica prezentarea considerând că toate semnalele au aceeaşi energie (ceeace nu e întotdeauna adevărat, decât dacă punctele ce reprezintă semnalele se află pe un cerc, cu centrul în origine) ; 0, 1, , -1iM==εε i Reamintim că apare o eroare, dacă statistica ataşată ipotezei adevărate, este depăşită de oricare altă statistică Probabilitatea de apariţie a unei astfel de erori se notează cu H {} {} 12 100max, , , MPPT TT T− =H >> > {} 12 100 0 0 MPT TsiT T si siT T− ∞ ==H >> > () {} ∫ , , , t ,MTPt T t T T T t p tdt−012 100 −∞ Deoarece statisticile sunt statistic independente, probabilitatea deciziei corecte devine ∞ , , , t ,P PtTtT T T t ptdt==>> >H () {} ∫ 012 100MT− −∞ 1M∞ − ⎡⎤ = ,Pt TT t p tdt=>H () {} ∏ ∫ 000Ti⎢⎥ 1i−∞ = ⎣⎦ în care 2 ⎫ 21t⎧ε − () ⎪⎪ exppt=− () ⎬ 022T⎨ ε 2σ ε 2πσ ⎪⎪ ⎩⎭ Se știe că, pentru o variabilă aleatoare X, repartizată normal, adică 2 Xμσ∼N, () se poate scrie că 2t ⎧⎫ − () 1xtμμ − ⎪⎪ ⎛⎞ − =Φ {} ⎨⎬ ⎜⎟ 2expPX t dx H {} {} 12 100max, , , MPPT TT T− 12M− ∞ ⎫ − ε 1221tt⎧ ⎡⎤ () + ε ⎪⎪ ⎛⎞ − ⎬ ⎜⎟ ⎢⎥ expdt=Φ2 ∫ 2⎨ 22σ ε ε σ ⎝⎠ ε πσ ⎣⎦ −∞ ⎪⎪ ⎩⎭ Dar probabilitatea de eroare adunată cu probabilitatea deciziei corecte dă unu, așa că = {} {} 1 ,12 100max , , MPT T T T− 12M− ∞ ⎫ 1221tt⎧ ε − ⎡⎤ () + ε ⎪⎪ ⎛⎞ Φ − ⎬ ⎜⎟ ⎢⎥ 1 expdt=−2 ∫ 2⎨ 2σ ε ε ⎝⎠ 2σ πσ ε ⎣⎦ −∞ ⎪⎪ ⎩⎭ 15 Cu schimbarea de variabilă 2t+ε u= ε σ expresia probabilității de eroare devine 12M− ∞ ⎛⎞ 1⎧⎫ε ⎪⎪ Φ − − ⎡⎤ () ⎜⎟⎨⎬ e21expPuudu=− ⎣⎦ ∫ 2σ −∞ ⎝⎠ ⎪⎪ ⎩⎭ 1M− 2∞ 1⎧ ⎫ Φ − −⎡⎤ 1expuuENRdu=− () ()⎨⎬ ∫ 2⎣⎦ ⎩⎭ −∞ La un număr dat de simboluri transmise, M, de obicei o putere a lui doi, probabilitatea de eroare este o funcție descrescătoare de raportul ENR, definit prin ε 2ENR= σ Probabilitatea de eroare scade odată cu creșterea numărului M de simboluri, dacă ENR este același (vezi figura), deoarece trebuie să distingem între semnale a căror distanță între ele nu crește 16 Pentru transmisii se recurge, așa cum am mai afirmat, la semnale ortogonale Este necesar să lucrăm cu NM≥ În figură se dau exemple de semnale ortogonale, în cazurile N=M=2 și N=M=3, ce au energia unitară În cazul general, pătratul distanței (normei) de la punctul de semnal la origine, este energia semnalului Pentru cazul semnalelor de energie egală, punctele lor reprezentative trebuie să se afle pe un cerc, pe o sferă sau o hipersferă Cazul modelului liniar În cadrul modelului liniar vectorul de date, x, poate fi pus sub forma =+=+xHθ wsw w 0C∼N; , () Distincţia între două ipoteze se face prin valorile parametrului vector θ θ0H 0: ; semnalul util este absent= θθH 11: ; semnalul util este prezent= Problema de detecţie constă în a decide dacă semnalul util s este prezent în datele x Cele două ipoteze pot fi reformulate xwH 0: ; semnalul util este absent= xHθwH 11: ; semnalul util este prezent=+ Detectorul NP ia decizii conform regulii 1H 11TT−− ′ ==>xxCsxCHθ () 1Tγ Asigurând, la o probabilitate impusă a alarmei false, o probabilitate de detecţie dată de 11T−− sC s () ) ( DFAPQQP=− 11TT−− θ HC Hθ () ) ( FAQQ P=−11 17 Exempu Detecţia unei sinusoide din zgomot alb, gaussian Semnalul util este de forma φφπ φπ=+= − 000[ ] cos(2 ) cos cos 2 sin sin 2sn A fn A fn A fnπ ab =+ ∈ () cos 2 sin 2 ;000 0, 0 5afnbfnfππ Cele două ipoteze se pot formula astfel H 0: [ ] [ ]; 0, 1, , -1xn wn n N== =++=H 100: [ ] cos 2 sin 2 [ ]; 0, 1, , -1xna fnb fnwnn Nππ 2 wnσ [ ] 0, ∼N () Pentru ipoteza H1 putem scrie că 10 xw⎡⎤ ⎡⎤ ⎡ ⎤ ⎢⎥ ⎢⎥ ⎢ ⎥ ⎡⎤ 00cos 2 sin 2 ffxawππ ⎢⎥ ⎢⎥ ⎢ ⎥ =+ ⎢⎥ b⎢⎥ ⎢⎥ ⎢ ⎥ ⎣⎦ ⎢⎥ ⎢⎥ ⎢ ⎥ θ 1cos 1 2 sin 1 2 NfNfxN wNππ −− () () 00−− ⎣⎦ ⎣ ⎦ ⎣⎦ xwH xHθw =+1 Pentru zgomotul alb gaussian, matricea de covarianţă şi inversa ei sunt 121− ==CI C I uuσ 2 σ Decidem că semnalul util este prezent, dacă 1H 11TT− ′ ==>xxCHθxHθ () 112Tγ σ Statistica ce se utilizează în luarea deciziei se poate transforma 1H 2 22σ T ′′ == >xxxHθ ()() 1TTγ NN Expresia acestei noi statistici este T 22⎛⎞ T T == HHxxxθθ () 11T′ ⎜⎟ NN⎝⎠ Expresia marcată cu roşu este estimatorul parametrului θ1 18 x⎡ ⎤ ⎢⎥ − ⎡⎤ () 001 cos 2 cos 1 2 22TfNfxππ ⎢⎥ Hx = ⎢⎥ − ⎢⎥ () 000 sin 2 sin 1 2fNfNNππ ⎣⎦ ⎢⎥ xN− ⎣⎦ 1N− 2⎡⎤ ∑ 0[]cos2xn fnπ ˆa⎢⎥ N= 0n⎡⎤ ˆ⎢⎥ θ = ⎢⎥ N==1 1− ˆ 2b⎢⎥ ⎣⎦ ∑ 0[]sin2xn fnπ ⎢⎥ 0nN= ⎣⎦ Noua statistică devine a⎡⎤ ˆˆˆT ⎡⎤ ==+=xθθ () 11ˆˆTabaabb′ ⎢⎥ b⎣⎦ ⎣⎦ şi este o corelaţie între vectorul estimat din date, cu o replică păstrată la recepţie *** Pentru a generaliza rezultatul obţinut în exemplul prezentat, plecăm de la expresia estimatorului MVU pentru modelul liniar 1− 11TT− ˆ− =θHC H HC x () şi punem statistica detectorului NP pentru modelul liniar, sub forma 11TTT−− xxCHθHC x θ () () 11T== Modificăm această formă a statisticii, prin introducerea a doi factori cu produsul unitar 11TTT−− xxCHθHC x θ () () 11T== T ⎡⎤ 1− −− 111TTT− ⎢⎥ HHC H HCx θCH ()() =1 ⎢⎥ θ 1ˆ⎢⎥ ⎣⎦ şi rezultă 1TT− ˆ x θ HC Hθ () 11T= Matricea de covarinţă a estimatorului, şi inversa ei sunt 1− 111TT− −− CHCH CHCH () ˆˆ;== θθ 19 Statistica testului ia forma finală 1− ˆ x θC θ () ˆ11T= θ Deciziile se iau conform regulii modificate 1H 1− ˆ ′ =>x θC θ () ˆ11Tγ θ Dacă analizăm comparativ statistica obţinută acum, cu cea din relaţia (14) H 11H 11T−− ˆ ′ γ==>>x θCxθxCs ()() ′ TTγˆ11;; θ rezultă că se pot menţine cele stabilite pentru un semnal de replică cunoscut şi în cazul parametrului de replică cunoscut, făcând substituirile formale ˆ →θCxsθC , ˆ1,→→ θ Problema de detecţie se transformă H θθ0s0sH == ⎫ ⎧ 0000::;;== → ⎬⎨ H sθsθ H = ⎩ =1111:;:; ⎭ Probabilitatea de eroare este, tot ca urmare a celor remarcate 1111T−−−− =→−−θCsCθs ()() ˆ11DADFAFPQQPPPQQ= ) ) ( ( θ Vom da un rezultat important, fără a-l justifica: Pentru a testa ipoteze despre un parametruθ al unui model liniar, vom înlocui cele N eşantioane de date cu p valori estimate cu 1− 11TT− ˆ− =θHC H HC x () Repartiţia vectorului estimator este 1− 1T− ˆ HC H∼N , ;θθ () () Orice problemă de decizie privind vectorul θ, se poate rezolva apelând la estimatorul dat 20 21